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Abstract: There are different schemes based on observers to detect and isolate faults in 
dynamic processes. In the case of fault diagnosis in instruments (FDI) there are different 
diagnosis schemes based on the number of observers: the Simplified Observer Scheme 
(SOS) only requires one observer, uses all the inputs and only one output, detecting faults 
in one detector; the Dedicated Observer Scheme (DOS), which again uses all the inputs 
and just one output, but this time there is a bank of observers capable of locating multiple 
faults in sensors, and the Generalized Observer Scheme (GOS) which involves a reduced 
bank of observers, where each observer uses all the inputs and m-1 outputs, and allows the 
localization of unique faults. This work proposes a new scheme named Simplified Interval 
Observer SIOS-FDI, which does not requires the measurement of any input and just with 
just one output allows the detection of unique faults in sensors and because it does not 
require any input, it simplifies in an important way the diagnosis of faults in processes in 
which it is difficult to measure all the inputs, as in the case of biologic reactors. 
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1. Introduction 

Processes supervision systems for operators have evolved as new techniques of detection and 
isolation of faults have appeared. Research in this field has also grown as the complexity of industrial 
processes has grown and this has motivated the development of different focuses for FDI 
system design. 

The diagnosis of faults can be done using observers. One great advantage of the diagnosis schemes 
based on observers is that in comparison with other methods they are very large schemes. The high 
level of complexity in current industrial processes has led to a situation where the amount of 
information generated by these processes can overcome the capacity of analysis of human operators, 
which hinders decision making [1,2]. The most recent supervision systems have the capacity to carry 
out diagnosis and maintenance functions, in order to guarantee the correct functionality of highly 
complex processes [3-8]. 

As an example of a process that is hard to supervise, in this work the production of biogas in an 
anaerobic reactor is used as a case study in which faults are diagnosed and isolated using a scheme 
based on observers. In many publications about non linear observers for the design of FDI systems, the 
residuals are based in the error of the estimation obtained by the observer [9]. In biological processes, 
due to their non linear nature, in the majority of cases they are not completely observable, therefore it 
is more appropriate to consider some relationships among parameters, instead of attempting to estimate 
them individually [10,1 1]. The work presented in [1 1] explores a methodology to determine the global 
state and parameters of biological reactors. The method proposed in the article allows one to formalize 
the design of asymptotic observers, which are capable of evaluating certain variables of state, which 
are not measured for the anaerobic digestion process, despite certain doubts about the kinetics of the 
process. We should mention that asymptotic observers need the information of all the input variables 
of the process, which it is practically impossible to achieve under normal conditions in the operation of 
an anaerobic reactor. Overcoming this problem implies knowing all the input variables of an anaerobic 
reactor and for this purpose an observer by intervals was developed in [12]. The main characteristic of 
intervals observers is that they are capable of providing estimations of guaranteed intervals of non 
measured state variables instead of an exact estimation, if a superior and inferior limit is provided for 
each one of the input processes. 

2. Fault Diagnosis Schemes 

Diagnosis schemes based on observers can be classified according the type of fault detected: sensor 
faults (Instrument Fault Detection or IFD), actuator faults (Actuator Fault Detection or AFD), and 
component faults (Component Fault Detection or CFD). Diagnosis schemes can also be classified 
according the number of observers that are used. There are schemes with one observer: a Direct 
scheme is a scheme of just one observer of complete order. The Simplified Observer Scheme (SOS), is 
a scheme of one observer of reduced order. For sensor faults (IFD), the only observer in this scheme 
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uses all the inputs and one output, which only provides simple redundancy and only allows the 
localization of faults in one sensor. In AFD, the only observer uses all the outputs and just one input. 
When several observers constitute a bank of observers of reduced order we have a Dedicated Observer 
Scheme (DOS). For faults in sensors (IFD), each observer uses all the inputs and just one output. The 
number of observers equals the number of outputs (sensors). For actuator faults (AFD) each observer 
uses one input and all the outputs. It should be mentioned that the DOS scheme allows the localization 
of multiple faults, either in sensors (IFD) or in actuators (AFD). The Generalized Observer Scheme 
(GOS) is formed by a bank of observers of reduced order. For faults in sensors (IFD), each observer 
uses all the inputs and m-1 outputs, where m is the number of outputs. For actuator faults (AFD), each 
observer uses all the outputs and n-1 inputs, n being the number of inputs. 

3. Design and Implementation of a New Diagnosis Scheme: SIOS-IFD 

The SIOS-IFD is a scheme with just one interval observator, of reduced order, for faults in sensors. 
The main advantage of the SIOS-IFD scheme over all the previously presented schemes, is the fact 
that no input measurements are is required; it is only necessary to have knowledge of the interval of 
values that the named inputs can reach. SIOS-IFD only allows the localization of faults in one sensor, 
because it requires the in line measurement of just one output. 

Figure 1 shows a block diagram of the SIOS-IFD. In this figure it can be observed that the 
SIOS-IFD does not use the inputs Ui (I = 1,2,3,4, ...,n), but rather just uses one output (ys) to estimate 
the other two outputs {y 1 e y 2 ) an d in that way be able to generate the responses: r 1 =y 1 —$ 1 y 

r 2 = y 2 -y 2 - 



4. Mathematical Model of the Anaerobic Reactor Used as Case Study 

Next the simplified ADM1 mathematic model of the UASB reactor of the Institute Tecnologico de 
Orizaba, Veracruz, Mexico (Figure 2), which is the case study of the project of this article, 
is presented. 



Figure 1. SIOS-IFD. 
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Figure 2. UASB reactor. 




Simplified ADM 1 Model: 

s i 

x ± = Y 1 K m>1 — ——I pH x 1 — aD(t)x 1 — K d x 1 
A s,l s i 

s\ = D(s 1 — 5 X ) — K ml — ■ IpH x l 

QCH4 = (1 — Yi)Y CH4 K ml - 7— IpH^i — QcH4 (1) 

with: 

1 + 2 * 10°- s (p h ll-pHul) 

l V H = 1 + 10(PH- P H UL ) + IQ(PH LL - P H) 

where x\ is the concentration of the anaerobic mass, s\ is the concentration of organic matter expressed 
as chemic oxygen demand (COD), Qch 4 is the exiting flux of methane biogas, D(t) is the rate of 
dilution, K m i, Kd y K s j are the specific rates of growth of anaerobic mass, the dilution rate of the 
anaerobic reactor and the constant decrease of semisaturation for the anaerobic biomass, respectively. 
Yj is the coefficient of performance for the degradation of COD, s[ is the concentration of COD in the 
affluent, Iph represents the pH inhibition, where pH LL and pH UL are the lower and higher pH limits, 
respectively. The values of the model parameters are shown in Table 1 . 



Table 1. Model Parameters. 



Parameter 


Value and units 


K m l 


5.1 gCOD/gCODd 


K s l 


0.5 gCOD/L 


K d 


0.02 L/d 


Yj 


0.1 gCOD/gCOD 


YcH4 


0.35 LCH 4 /gCOD 


a 


0.5 (adimensional) 
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5. Interval Observer Designed for the SIOS-IFD Scheme 



In this section the design of the interval observer designed for the IOS - IFD scheme is presented. 
The designed interval observer is capable of stimating value xi y Qch 4 from the in-line measurement of 
si.lt should be emphasized that the designed interval observer does not require any measurement of 
the reactor input variables. The first necessary condition for the design of an interval observer is that a 
hypothetic observer of known inputs must exist, called base observer. In order to satisfy this first 
condition, the asymptotic observer presented next was designed. 

The model described by the set of differential non linear Equation (1) can be rewritten in the 
following way: 

m = [*i Qch 4 StY = Cf(x(t),t)+A(t)x(t) + b(t) (2) 

with: 

/(x(t),t) = [A*i*i] 

C = [Cl \ CU T = [Yi a~Y 1 )Y CH4 I -If 

[-(aD(t) + k d ) 0 

A(t) = 



A xl {t) 
A 2 i(t) 



0 



0 



-1 
0 



0 
0 

-D(t). 



b(t) = [bl(t) \ bUt)]=D(t)[0 0 isKOl 1 



where x{f) E 9T is the state vector, C E$i mXn is the matrix of performance coefficients and 
f(x(7), t) E yi m is the vector that contains the non linearity of the model, which are assumed to be totally 
unknown, the time variant matrix A(t) E $l nXn is the matrix of state and b(t) E 9T is the vector of the 
observer entries. 

The asymptotic observer is designed under the assumption that all inputs are known, and m 
measures states on-line. Thereby the the space of states can be divided in such way that Equation (2) 
can be rewritten as: 

7 (0 = Cjf(x(t),t) +A 11 (t)v 1 (t) +A 12 (t)v 2 (t) + hit) 

2 {t) = C 2 f(x(t),t) + A 2 i(t)v 1 (t) + A 22 (i)v 2 (i) + b 2 (t) 

where the m measured states are regrouped in vector v 2 (t) (dim v 2 (t) = m) and the variables that will be 
estimated are represented by vi(t) (dim vj(t) = s = n - m). The matrices An(t) E $l sX \ An{t) E 9l* Xm , 
A lx ii) E 9T X *, ^22(0 E 9T Xm , Ci E 9T r , C 2 E 9T Xr , b x E W y b 2 E W n are the corresponding partitions of 
A(t), Cy b(t) 9 respectively. 

Equation (3) represents the asymptonic observer that was designed for the anaerobic reactor 
described for the Equation (1): 

ft = W(t)w(t) + g>(t)v 2 (t) + Nb(i) 

w(0) = NMO) + N 2 v 2 (0) (3) 
x 1 (t) = N^\w<it)-N 2 v 2 (t)) 
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with: 

w(t) = (JV^nCt) + iM 21 (t))JVf 1 

*(t) = N t A 12 (f) + N 2 A 22 (t) - W(t)N 2 

In the design of this asymptonic observer, it was assumed that the rate of dilution D, and the inputs 
to the digestor are considered as knowntherefore A(t) y b(t) are known V t > 0. The asymptonic 
observer 1 was designed assuming si as the only state measured on-line, meaning, V2(t) = [si] (dim 
V2(t) = 1), while xjy si were assumed as estimated states, meaning vj(t) = [xj Qch 4 ] T (dim vj(t) = 2). 
For the calculation of the matrix N, Nj was arbitrarily chosen as an identity matrix, to calculate N2 in 
the following way: Ni= -N1C1C2K (where C2 § is the pseudo reverse generalization of C2), thus 
obtaining the following result: 

10; Y x 
0 1 ! (l-Y t )Y CH 4\ 

Substituting Ni and N2 in Equation (3), the following is obtained: 



n = [n 1 = yv 2 ] = 



W (t) = w e = [-( aD V + k -) _°J 



Finally, W e + was calculated using the minimum value of the dilution rate D(t) = 0.01, obtaining: 



r + _ r-0.025 0 1 

e ~ y 0 -iJ 



W e ' was calculated using the maximun value of the dilution value D(t) = 1, obtaining: 



e 



= [-0.52 0 1 
I 0 -V 



The eigenvalues of W e + y W e ~ are, respectively: 

W> [-0.025 -1] 

Mw;) = [-o.52 -i] 

Now, in order to make the observer asymptotically stable, the following conditions must be 
fulfilled: 

(a) W eiJ >0Vi^j. 

(b) We y W e ~ must be stable. 

Since both conditions are fulfilled, the observer is asymptotically stable. Next the design of an 
interval observer based on the designed asymptotic observer is presented. Once the asymptotic 
observer that will work as base observer have been designed, we continued with the design of the 
interval observer, which is based in the supposition that the values of the input vector b(t) are 
unknown, but their upper limitb + (t) and lower limit b~(t) are known, so b~(t) < b(t) < b + (t). 
Thus, the asymptotic observer presented in Equation (3) works as a structure base to build the next 
interval observer: 
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For the upper limit: 



(ti + (t) = Ww + (t) + 0(t)v 2 (t) + M£ + (t) 
w + (0) = 7V 1 v 1 + (0) + 7V 2 v 2 (0) 
x 1 \t) = N 1 -\w + (t)-N 2 v 2 (t)) 



(4) 



Ly 2 + (0. 



For the lower limit: 



w-(0) = 7V 1 v 1 -(0)+7V 2 v 2 (0) 
xr(t) = N 1 - 1 (^"(t)-Af 2 i; 2 (t)) 



(5) 



y 2 (0 



where: 
M = [N t 

P + (t) = 



n 2 i 



Kit) i -(fc 2 + (t) + & 2 (t)) 



*r(0 i 2 (&2+(t) + 62_(t)) 



-(& 2 + (t)-& 2 -(t)) 



-2^2 (0-^(0) 



The convergence of the interval observer is based in the principal of cooperation defined by [13] 
(see Lema 1). Let e + (t) = v t + — v x y e~(t) =v 1 — v{~ be the errors of estimation associated with 
Equations (4) and (5), respectively. The follwoing expression represents the dynamics of those errors: 

§*(t) = W e e*(t)+f2* (6) 

It uses to refer to any of the errors e + (t) o e~(t) so that their dynamics have the same 
mathematical structure. In the previous equation f2* = N 1 ~ 1 N( < b + (t) — b(t)) for the case of the upper 
limit, while £T = N 1 ~ 1 N{b(t) - b~(t)) for the case of lower limit 
Lema 1 [11] 

In a) = f(a), t). The system it is cooperative if ^L^—il > q V/i ^ d. 

da) d 

This implies that if o)(0) > 0 then > 0 Vt > 0. Therefore, if the initial conditions of the 
estimated variables are unknown but have their limits organized as: #i~(0) < v ± (0) < #i + (0) then 
§*(0) > 0 . Besides, if the system (6) is cooperative, the matrix W e is stable (Hurwitz) and 
if Q* is positive (or zero), it can be guaranteed that §*(t)>0Vt>0, and in consequence: 
0i~(t) < v ± (f) < v^ifiVt > 0 [12]. 



6. Experimental Results of the Interval Observer 



To experimentally verify the operation of the developed interval observer, the anaerobic digester 
was fed with wastewater from a brewery, which had maximum values of 3 gCOD/L and minimum 
values of 2 g COD/L, these values are based to the digester being accustomed to consume 3 gCOD/L. 
If the concentration of the water is greater than this, it can be diluted to achieve the desired value, but 



Sensors 2011, 11 



619 



if it is lesser, it will be difficult to remove excess water to achieve the desired value, being in the worst 
case a concentration of 2 gCOD/L. The dilution rate is bounded by a maximum value of 0.74 d 1 and 
minimum value of 0.26 d 1 (these parameters were obtained experimentally by applying to the plant a 
positive bounded control as explained by Zavala [14], which sets the value of D by reference to the 
values of 3 g COD/L and 2 g COD/L for SI, these being the maximum and minimum value that is 
generated from D). The concentration of organic matter in the effluent (si) was measured on-line to 
estimate xl and Qch 4 > Figure 3 shows the convergence of the range observer for the concentration of 
anaerobic mass xl; Figure 4 shows the curves, and on-line measurement for the methane biogas 
flow Qch 4 - 

Figure 3. Estimation of the concentration of anaerobic mass xl. 
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Figure 4. Estimation of methane biogas flow Qch 4 - 
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7. Experimental Results of the SIOS-IFD 



The developed SIOS-IFD scheme is able to detect unique, sudden and permanent faults in sensors 
si, xi and Qch 4 , from response r x = y 1 — 5^ and response r 2 = y 2 — y2> being yj = xj and y2 = Qch 4 - 
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Figure 5 shows the responses to a single, sudden and permanent fault, +5% in the sensor si, on day 
45 of experimentation. This figure shows that both responses have faults on day 45, because both 
x ± and Q C H4> have been estimated from si so the two estimates are wrong from the moment in which 
the fault of si occurs. 

Figure 5. Responses to a single sudden and permanent fault, +5% in the sensor sj. 



5 
4 
3 

I'X1 2 
1 

0 



1QCH4 



10 



20 



3D 
Days 



40 



60 



























































0 ] 


L0 2 


O 30 A 


0 50 60 



Days 



Figure 6 shows that each of the responses generated by the diagnostic scheme developed are filtered 
through a moving window averaging filter. Subsequently, each filteredresponse is evaluated by a 
threshold detector. The output of each threshold detector is connected to a display, which can only 
display one of the following values: -1, 0 or 1. The value "-1" indicates that the response reached a 
value averaged less than the threshold, so, the response had a negative average change. Similarly, the 
value "1" indicates that the response had an average positive change. Finally, the value "0" indicates 
that either the responses does not have any media change; or rather variations in their average values 
were below preset thresholds. 

The array of symptoms presented in Figure 6 [11] is the result of the evaluation of the responses 
generated by the fault diagnosis scheme developed for the example of a single, sudden and permanent 
fault of +5% in sensors;. These responses were shown in Figure 5. 

Figure 6. Developed response assessment scheme. 
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Table 2 presents the Structured Diagnostic Matrix for unique faults in the sensors si, x\ y Qch p 
obtained from the developed response assessment scheme. 

Table 2. Structured Diagnostic Matrix. 





Faults/ 


Faults/ 


Fault*/ 


Fault*/ 


Faultg«f 4 


Faulty 




+5% 


-5% 


+5% 


-5% 


+5% 


-5% 


ri 
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-1 
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-1 
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-1 
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Shaded In Table 2 is symptom vector [11], which is the result of the evaluation of the response 
generated by the fault diagnosis scheme developed for the example of a single, sudden and permanent 
fault of +5% in the sensor sj. In this same Table it can be seen that there is a different symptom vector 
for each single, sudden and permanent fault in sensors s^xiy Qch 4 - 

8. Conclusions 

In this article experimental and simulation results of a novel system of diagnosis of faults in sensors 
are presented. It has been named as SIOS-IFD (Simplified Observer Interval Scheme — Instrument 
Fault Detection). The main advantage of the SIOS-IFD scheme in comparison with all the other 
schemes presented above is the fact that the measurement does not require any input, only knowledge 
of the range of values that such targets can achieve is required. The SIOS-IFD scheme only allows the 
location of a single sensor fault, for which it requires the in-line measurement in a single output. 
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